Formation of Liesegang patterns: A spinodal decomposition scenario 
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Spinodal decomposition in the presence of a moving particle source is proposed as a mechanism 
for the formation of Liesegang bands. This mechanism yields a sequence of band positions x n that 
obeys the spacing law x n ~ Q(l + p) n . The dependence of the parameters p and Q on the initial 
concentration of the reagents is determined and we find that the functional form of p is in agreement 
with the experimentally observed Matalon-Packter law. 
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Pattern-forming chemical, physical and biological pro- 
cesses are common in nature and patterns often emerge 
in the wake of a moving front pL In particular, when an 
electrolyte A diffuses into a gel containing another elec- 
trolyte B, the eventual formation of a rhythmic pattern 
of precipitate by the moving chemical reaction front is 
known as the Liesegang phenomenon . The observed 
precipitation patterns usually consist of a set of bands or 
rings (depending on the geometry of the system) clearly 
separated in the direction perpendicular to the motion of 
the front. This phenomenon is believed to be responsible 
for many precipitation patterns such as e.g. the structure 
of agate rocks || . Although the Liesegang phenomenon 
has been studied for a century, the mechanisms respon- 
sible for these structures is still under discussion Q] . 

Most of the reproducible Liesegang patterns are char- 
acterizable by the following generic laws. First, the po- 
sition of the n-th band x n (measured from the initial 
interface of the reagents) is proportional to \ft^ where 
t n is the time elapsed till the appearance of the band. 
This so-called time law || is a direct consequence of the 
diffusive nature of the dynamics. Secondly, the positions 
x n of the bands usually form a geometric series (spacing 
law §): 

n large 

x n >Q0-+pT (i) 

where p > is called the spacing coefficient and Q is the 
amplitude of the spacing law. Finally, the width w n of 
the bands have been observed to increase with n and to 
obey the width law 0, w n ~ x n . 

Most of the detailed experimental observations con- 
cern the spacing law. It has been found that the spac- 
ing coefficient p is a nonuniversal quantity depending 
(among other parameters) on the experimentally control- 
lable concentrations ao and bo of the outer (A) and inner 
(B) electrolytes. This dependence is expressed by the 
Matalon-Packter law ||||: 

p = F(6 )+G(6o)- , (2) 
a 

where F and G are decreasing functions of their argu- 
ment b JlOfl . 



The task of theories is to explain the existence of pat- 
terns obeying the spacing law (|l|) and, on a more sophis- 
ticated level, to derive the Matalon-Packter law @ . The 
theoretical approaches proposed up to now can be di- 
vided into two categories. The first one contains the ion- 
product super saturation type theories where the 
outer and inner reactants turn directly into precipitate 
(A + B — > D) whenever their local concentration prod- 
uct is above a threshold q* . The second category con- 
tains the intermediate- compound theories J^|,[]^| which 
assume the existence of a species C (A + B — > C — > D). 
The C-s are produced in a reaction-diffusion front that 
moves through the system and the precipitation C — > D 
takes place only if the local concentration c reaches some 
threshold value c*. In the presence of D, the process con- 
tinues until the concentration of C drops below another 
threshold d*. The problem with these theories is that 
either they contain parameters such as e.g. d* that are 
difficult to control experimentally and not easy to grasp 
theoretically, or they describe a detailed mechanism that 
is too complicated to deduce quantitative consequences 
such as the Matalon-Packter law. Our goal with this Let- 
ter is to describe the formation of Liesegang patterns in 
terms of a model that contains the basic mechanism of 
phase separation that underlies the band formation but, 
at the same time, simple enough so that quantitative 
predictions are readily made. In particular we show the 
existence of the spacing law, obtain the Matalon-Packter 
law with estimates of F(bo) and G(bo) and, furthermore, 
the bo /ao dependence of the amplitude Q of the spacing 
law is also determined. 

All the theories discussed above can produce the spac- 
ing law but the comparison with the Matalon-Packter 
law suggest that the intermediate-compound theories are 
clearly preferable Q. Accordingly, we shall assume that 
the motion of the reaction-diffusion front and the dynam- 
ics of the C particles deposited by the front should be the 
main ingredients of a theory of Liesegang patterns. The 
dynamics of the A + B — > C front has been solved so its 
motion and the production of C-s are well known [ jl6| . 
Thus the new aspect of our theory is the model for the 
dynamics of C-s. 
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In our picture, the C-s are particles that can move only 
by diffusion due to the presence of the gel. There is an 
attractive interaction between the particles that induces 
aggregation at low enough temperatures and high enough 
densities. An easily understandable, spatially discretized 
version of the model is the following spin- 1/2 kinetic Ising 
model with competing spin-flip and spin-exchange dy- 
namics Empty and occupied lattice sites are asso- 
ciated with down and up spins, respectively. The ini- 
tial state is empty (all spins are down) and the moving 
front flips the down spins at a given rate (Glauber dy- 
namics The diffusion is described by a spin ex- 
change process (Kawasaki dynamics The rates of 
exchanges arc governed by a heat bath at temperature 
T with ferromagnetic couplings between the spins mod- 
eling the attraction among the C"s. The experimentally 
observed freezing of the emerging patterns implies that 
the corresponding dynamics takes place at a very low ef- 
fective temperature. If this model represents the pattern 
forming process correctly then one expects the emergence 
of bands of up and down spins in the wake of the mov- 
ing spin-flip front. In order to understand how the bands 
arise, let us consider the phase diagram of an Ising model 
depicted on Fig.l. 
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FIG. 1. Qualitative phase diagram for the Ising model. 
The magnetization m is related to the density c of C par- 
ticles through m = 2c — 1. The solid line is the coexistence 
curve and the dotted one is the spinodal line. Si is the initial 
state with m = — 1, ±m e are the equilibrium magnetizations 
at a given temperature T while ±m s are the magnetizations 
at the spinodal line. The front alone would leave behind a 
magnetization mo- 



One starts from the all-spins-down state (empty state, 
Si in Fig.l). Since an A+B — > C reaction front leaves be- 
hind a constant density cq of C-s ||lq| , the spin-flip front 
is chosen such that it produces a local magnetization 
too = 2co — 1. We would like the front to bring the sys- 
tem into the unstable regime, thus too > m s is assumed. 
As to is increasing from — 1 to toq in the front , the local 



state moves from Si towards So (Fig.l). With time, one 
crosses the coexistence line (to = m e ) and enters into 
the metastable phase 1 20 . Small clusters of up-spins nu- 
cleate at and aggregate behind the front. However, the 
nucleation is an activated process and so its characteris- 
tic time-scale T nuc i is large at low temperatures. If T nuc i 
is much larger than the time Tf ront needed by the front 
to put out the local magnetization too then the system 
reaches the unstable state i.e. crosses the spinodal line 
(to = m s ) (|lj. Once the spinodal line is crossed, the 
phase separation takes place on a short time-scale and a 
spin-up domain is rapidly formed at or behind the front, 
hence the formation of a Liesegang band. 

This band acts as a sink for the up-spins and, in the 
vicinity of the band, the local magnetization decreases 
and the front is no longer in the unstable region of the 
phase space. When the front moves far enough, the effect 
of the band as a sink diminishes. Thus the magnetization 
grows and the spinodal line is crossed again resulting in 
the formation of the next band. The repetition of this 
process should lead then to the Liesegang pattern. The 
new feature of the picture described above is the assump- 
tion that the state of the front is quasiperiodically driven 
into the unstable regime. 

The above microscopic picture can be described on 
the mesoscopic level as follows. The diffusive dynamics 
of the coarse grained magnetization to(x, t) is described 
by Model B of critical dynamics |^2| with the moving 
spin-flip front appearing as a time-dependent source term 
S(x,t): 



dtm(x, t) 



Xdl 



em(x, t) — r ym 3 (x 1 1) + ad^.m(x, t) 



+ S(x,t) 



(3) 



where A is a kinetic coefficient, e measures the devia- 
tion from the critical temperature T c and e > ensures 
that T < T c . The parameter 7 is positive to guar- 
antee overall stability and a > provides the stabil- 
ity against short-wavelength fluctuations. In principle, 
equation (^) should contain two noise terms. First, there 
should be thermal noise that, in the absence of the source, 
would bring the system to equilibrium described by the 
phase diagram in Fig.l. Second, there should be noise in 
the source S since the origin of this term is a reaction- 
diffusion process. We shall omit both noise terms. The 
thermal noise is neglected since the effective tempera- 
ture is expected to be low as discussed above. As to the 
source term, it has been shown that the properties of the 
A + B — > C type reaction fronts are mean-field like above 
dimension two p3|| . We take this as an indication that 
the noise in S can be neglected. 

In order to complete the description of the model, we 
need to discuss the actual form of the source term, S. The 
problem of the A + B — > C reaction-diffusion front is well 
understood at the mean- field level fu|. The source term 
is the production rate of the C particles that can be cal- 
culated analytically. Provided the reagents A and B are 
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spatially separated at the initial moment [i.e. their densi- 
ties are given by a(x < 0, t = 0) = do, b(x < 0, t = 0) = 
and a(x > 0, t — 0) = 0, b(x > 0, t = 0) = b ], the source 
term is a gaussian to an excellent accuracy: 

s(xt) - A exp r fr-*/(*)] 2 i f4) 



The center of the front moves as ar/(£) = ^/2Dji with 
the diffusion constant D / given by the following equation 
eif(y/D f /2D) = (a - &o)/(a + bo) with D = D a = D b 
being the diffusion coefficient of the A and B particles 

The front is well localized since its width increases with 
a small power of time w(t) = 2y/Dt 1/6 /(fca if) 1/3 where 
k is the reaction rate of the A + B — > C process, and 
if is given by K = (1 + VaoHV^T 1 exp (-D f /D). 
Finally, the amplitude of the source can be expressed as 
A = 0.3fca§X 4 / 3 . Thus the initial densities (ao,&o), the 
diffusion constants (£> = D a = Df,), and the reaction 
rate k determine all the parameters in S(x,t) |l6| , |24| ]. 

We turn now to the solution of eqn (|3|). In the ab- 
sence of the source term, S = Q, the globally stable 
homogeneous solutions are m" = ±-\/ e /7- One can 
also find other homogeneous solutions m(x, t) = m^ 
which, however, are unstable to small perturbations if 
\ m h\ < vWI^t) = m s- The value m s gives the location 
of the spinodal line ||0[ . 
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FIG. 2. Magnetization profile obtained for the following 
values of the front parameters: Df — 21.72, wo = 4.54, and 
A = 0.181 with length, time, and the magnetization (con- 
centration) measured in units of \fae, a/(Xe 2 ), and \J e/7, 
respectively. The dotted line denotes the rate of local mag- 
netization increase due to the source, S, measured in units of 
Ae 5/2 /(7 1/2 o-) and magnified by a factor 2 ■ 10 5 . The dashed 
line is the magnetization at the spinodal line, m s = — 1/V3. 

In the presence of the source term S, equation (|^) 
can be solved numerically using the initial condition 
m(x, t) — —ml and starting the source at the origin. 
The source alone would leave behind a uniform particle 



density cq w 0.85!fao y/D/Df [|l6| corresponding to a 
magnetization mo = 2cq — 1. The solution that evolves 
depends crucially on the value of mo- 

If I mo I > m s , the system is always outside of the unsta- 
ble domain and a uniform band of precipitate is formed. 
Since there is no noise in the system, this band is sta- 
ble. Phase separation does take place when |mo| < m s 
and a pattern similar to that shown in Fig. 2 is observed. 
If mo is near the spinodal value then, at early stages of 
evolution, a nearly periodic set of narrow bands emerges 
that coarsens with time. Later on, the newly formed 
pattern crosses over to Liesegang type bands with the 
distance between consecutive bands increasing. If mo is 
such that the systems is deep in the unstable domain, 
well defined Liesegang bands form from the very begin- 
ning (Fig. 2 shows an example of this type of evolution). 
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FIG. 3. Positions of the bands a function of their 

order of appearance n. Length is measured in units of y/ae. 
The full line is the best fit given by equation (|H|). 

In agreement with experiments, once formed, the 
Liesegang type bands are static on the time-scale we are 
able to observe (occasionally up to 100 bands were gener- 
ated). Thus their spacing is well defined and can be stud- 
ied. Fig. 3 displays the positions of the bands of the pat- 
tern exhibited in Fig. 2. As one can see, a two-parameter 
fit of the form 

x n = Q[exp(pn) - 1] (5) 

gives excellent agreement with the data (similar qual- 
ity fits can be produced for a wide range of parame- 
ter values). Thus we see that the model contains the 
spacing law ([!]) with the spacing coefficient given by 
1 +p = exp(p). 

Next, one can investigate whether the patterns in our 
model obey the Matalon-Packter law. Calculating the 
p-s for a set of a at a fixed b , we find (see Fig. 4) that 
the results indeed agree with the linear 60 /&o dependence 
given by eqn.(||). The slopes and the intersections with 
the p axis of the straight line fits provide us the functions 
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F(bo) and G(bo). Due to the restricted interval of 60 val- 
ues that are available we can deduce only an approximate 
functional form for these functions. It is clear, that both 
F and G are decreasing functions of their arguments, in 
agreement with experimental findings ||[|. A fit to a 
power-law form gives F(6 ) ~ &o anc ^ G(bo) ~ b, 
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From the point of view of applications, it is important 
to know the ao and bo dependence not only for p but also 
for the amplitude Q of the spacing law. In our model 
we can investigate this quantity as well and our prelim- 
inary results indicate that Q(ao,6o) ~ (ao/^o) ' 4 - The 
increase of Q with ao is a somewhat surprising result. It 
means that, as far as the number of bands in a given in- 
terval is concerned, there is a competition between p and 
Q. When ao is increased, this number increases since 
p ~ O.Q 1 but it decreases due to the change in the ampli- 
tude Q ~ af] a . This means that in experiments where n 
is finite (typically n ~ 20—30) one could observe thinning 
of the bands as a result of the increase of the concentra- 
tion of the outer electrolyte. 
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FIG. 4. Spacing coefficient p as a function of the ratio 
bo/ao for several values of 60. The unit of bo is fixed by the 
considering the front parameters given in Fig. 2 and setting 
bo/ao = 0.01 (filled circle in this figure). 



In conclusion we have proposed a new scenario for the 
formation of Liesegang patterns based on a spinodal de- 
composition mechanism. Our approach has the advan- 
tage of involving only a small number of parameters and 
there is no need to introduce artificial thresholds. Our 
model yields the Matalon-Packter law and allows the cal- 
culation of both the spacing coefficient p(ao,bo) and the 
amplitude Q(ao, bo) of the spacing law. 

This simple scenario can be improved by studying the 
role of the noise either by adding it to the continuum 
model (^) or by carrying out Monte-Carlo simulations on 
the kinetic Ising model with competing dynamics as de- 
scribed above. Work along these lines is in progress p5|. 
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